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ABSTRACT 

Some fraction of the material ejected in a core collapse supernova explosion may remain bound 
to the compact remnant, and eventually turn around and fall back. We show that the late time 
(> days) power associated with the accretion of this "fallback" material may significantly affect the 
optical light curve, in some cases producing super-luminous or otherwise peculiar supernovae. We use 
spherically symmetric hydrodynamical models to estimate the accretion rate at late times for a range 
of progenitor masses and radii and explosion energies. The accretion rate onto the proto-neutron 
star or black hole decreases as M esc t~^/^ at late times, but its normalization can be significantly 
enhanced at low explosion energies, in very massive stars, or if a strong reverse shock wave forms at 
the helium/hydrogen interface in the progenitor. If the resulting super-Eddington accretion drives an 
outflow which thermalizes in the outgoing ejecta, the supernova debris will be re-energized at a time 
when photons can diffuse out efficiently. The resulting light curves are different and more diverse 
than previous fallback supernova models which ignored the input of accretion power and produced 
short-lived, dim transients. The possible outcomes when fallback accretion power is significant include 
super-luminous (> 10'*'* ergs s~*) Type II events of both short and long durations, as well as luminous 
Type I events from compact stars that may have experienced significant mass loss. Accretion power 
may unbind the remaining infalling material, causing a sudden decrease in the brightness of some long 
duration Type II events. This scenario may be relevant for explaining some of the recently discovered 
classes of peculiar and rare supernovae. 

Subject headings: supernovae: general — supernovae: individual (SN 2008es, SN 1998bw, SN 2010X) 
— stars: massive — stars: neutron — accretion, accretion disks — black hole 
physics 



1. INTRODUCTION 

Ongoing optical surveys have discovered new classes 
of supernovae (SNe), incl uding sub-luminous (e.g., 
iPerets et al.ll2009l : lFolev et a l. 2009; Kasliwal ct al. 201^ 
and s uper-luminous (e.g.. iQuimbv et al.ll2011l : iGal-YamI 
|2012| ) events. Many of these events are difficult to 
explain in the context of standard models for which 
radioactive decay powers the optical light curve (e.g., 
SN 2008es, iGezari et ahl I2OO90 . Ahernative explana- 
tions fall into two main categories. Overluminous 
Type Iln SNe with narrow hydrogen absorption lines are 
thought to be powered, at least in part, by the inter- 
action of the supernova ejecta with circumstellar ma- 
terial, effectively re-thermalizing the supernova shock 
enerev (IS^th fc McCra vl[2007l: IChevalier fc Irwinl l2QlTI: 
iMoriva et all 12012: Ga l-YamI I2012D . Alternatively, the 
SN ejecta may be re-energized by the spindown power 
of a rapidly r otating magnetar w hich formed in the core 
collapse (Kas en fc Bildstenll2010l hereafter KBIO). For 
either of these two mechanisms to significantly modify 
the SN light curve, the energy input must occur at rel- 
atively late times (weeks to months) when radiative dif- 
fusion through the ejecta is efficient. 

Accretion onto a central compact remnant repre- 
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sents another potential means of injecting large amounts 
of energy in either successful or failed supernova ex- 
plosions. Compact object accretion i s associated 
with large-scale outflows in neutron stars (jFender et al.l 
20041) and stellar mass bla ck holes (microquasars, e.g., 
Mirabel fc Rodriguezll998D , and these outflows can carry 
away as much as ~ 10% of the gravitational bind- 
ing energy of the infalling material. This is particu- 
larly true when the accretion flow is rotationally sup- 
ported and radiatively ineffi cient ()Naravan fc Yil [1994 : 
IBlandford fc Begelmanl[T999l) . 

In "failed" SNe, the entire star accretes onto the cen- 
tral remnant, a black hole. If the progenitor lacks suf- 
ficient angular momentum to form a disk and hence 
tap the available accretion ene rgy, these events are "un- 
novae" (|Kochanek et al.l [2008D . i.e., stars disappearing 
suddenly from view. In the opposite case where even 
the iron core becomes rotationally supported, the ac- 
cretion energy may power a long-duration gamma ray 
burst (GRB) (the coUapsar m echanism, Wooslev 1993; 
iMacFadven fc WooslevI Il999( ). Much Ion ger gamma 
ray transients may also be possible if either the man- 
tle (IMacFadven et al.l 1200 ID or the hydrog en envelope 
(jWooslev fc Hegeiil2011t lOuataert fc Kasenil20T2l here- 
after QK12) becomes rotationally supported and drives a 
relativistic jet. The timescale associated with the energy 
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injection corresponds roughly to the free-fall time of a 
stellar layer, about ~ 0.1 s for the iron core, but as long as 
'-^ 1 yr for the hydrogen envelope of a red giant. Powerful 
winds from the accretion disk may eventually provide suf- 
ficient energy to turn the failed SN into a succ essful one, 
exploding the remainder of the star (Milosa vlievic et al.l 
I2OIOI: ILindner et al.ll2QH . 

In successful SNe, accretion from the "fallback" of the 
fraction of material remaing bound can be significant 
as well. The fallback may influence the resulting nu- 
cleosynthesis (Colgate 1971) or delay the pulsar mecha- 
nism in a you ng proto- neutron star such as in SN 1987A 
(|Michell [19881) . Early time fallback may also provide 
a link between the expl osion mechanism and the rem- 
nant mass distribution (iFrver et aLl I2012D or alter the 
radiated neutrino spectrum (jFrved 120091 ) . For red su- 
pergiant (RSG) progenitors with typical explosion ener- 
gies (~ lO^-'^erg), the fallback mass is relatively small 
(~ 0.1 Mq). However, in more compact stars (e.g., 
blue supergiants, BSGs) the formation of a strong reverse 
shock at the H/He interface can decelera te the ejecta and 
enhance the fallback mass to ~ 2Mq (IChevaliei<ll98l 
iZhang et al.l 1200^ . For weak explosions, most of the 
star may fall back, w ith only a small frac tion of the mass 
ejected in a dim SN (|Moriva et al.ll2010f ). 

While the dynamics of supernova fallback have been 
studied in a various contexts, little has been said about 
how fallback may impact the optical SN light curve. The 
energy released from fallback accretion may profoundly 
affect what we observe, if two conditions are met. First, 
the accretion energy must be injected at relatively late 
times (> days) otherwise it will be largely degraded by 
adiabatic expansion. Such late time accretion may be 
possible for progenitors with extended envelopes, or for 
those where a reverse shock develops and gradually slows 
the inner layers of ejecta. Second, the accretion energy 
must be thermalized within the SN ejecta. This is likely 
to occur if the energy injection takes the form of a nearly 
isotropic disk wind. If, on the other hand, the energy is 
in a beamed relativistic jet, we must consider whether 
the jet can breakout of the ejecta (perhaps producing a 
GRB) or whether it is trapped and thermalized in the 
interior. When these two conditions are met, fallback 
should produce a peculiar optical light curve, powered 
directly by the accretion energy. 

We study the impact of late time fallback accretion of 
SN light curves, and suggest that the wide range of po- 
tential events - from sub luminous to super-luminous - 
may be of relevance in explaining recent observations of 
peculiar SNe. In ^J2l we crudely estimate the efficiency of 
fallback-accretion-driven outflows. We numerically cal- 
culate accretion rates for a wide range of stellar pro- 
genitors to explore the variety of outcomes for accretion 
powered supernova light curves O, including sample 
light curves and their comparison to some recent unusual 
events f M.lj) . We also attempt to address the various 
requirements for these events to occur in Nature: the in- 
teraction of the outflows with both the infalling material 
and the outgoing ejecta fi j5.ll) . and angular momentum 
and disk formation ( ii5.2p . The major results are sum- 
marized in fj6l 

2. ACCRETION ENERGY 



At both low and high accretion rates compared to 
Lcdd/c^, where Lcdd is the Eddington luminosity, ac- 
cretion flows onto compact objects become hot and 
geometrically thick due to their inability to cool effi- 
ciently (icooi > ^infaii)- Such radiativcly inefficient ac- 
cretio n flows are expected to produce l arge-scale out- 
flows (INaravan k Yilll99l iBl andford & Begelma,nl 



Igumenshche v & Abra mowiczi i2000i: ,Pen et al.l 



1999; 



2003 



Begelmanl l2012: McKin nev et al.ll2012l) and/ o r Poyntini 



m nev et ai.llzUl^ and/ o r royntmg 
flux dominated jets (|De Villiers et al.1 120051 : iMcKinnevI 
|2006( ). This behavior is observed in the accretion flow 
onto the Galactic center black hole, Sagittarius A*, 
where the accretion rate at the Bondi radius (e.g., 
iBaganoff et al.|[2003t lQuataerj|2004[ ) is several orders of 
magnitude larger th an that onto the black hole (e.g., 
iMarrone et al.l [20071 ) . The fallback accretion rate fol- 
lowing a successful supernova explosion is highly super- 
Eddington and extremely optically thick to photons. For 
all timescales of interest here (> 1000s after the explo- 
sion) , the disk is not dense enough to cool by neutrino 
emission ([Kohri et al.l 120 05). We expect then that it 
should be radiatively inefficient, geometrically thick, and 
should drive large-scale outflows. 

The resulting mass out flow rate can be estimated fol- 
lowing [KDhrT^Fal] (|2005D by assuming that the accretion 
rate increases as some power of radius 



M(r) 



V^'fb 



(1) 



where Mfb and rfb are the accretion rate and radius at 
the outer disk edge, and < s < 1. We will write the 
radius in units of the Schwarzchild radius, r = Rc^/GM. 
The outflow speed should be comparable to the escape 
speed, Viu ~ c/ V2r, and the energy in each disk annulus 
is. 



dE^=^dM- = --^dr, 



(2) 



where ^ parameterizes our ignorance of the outflow 
physics, such as the fraction of fallback mass that is 
blown out again. The actual value of s is highly 



uncertain, but s 



1 is a reasonable choice 



, (e.g., 

Igumenshchev fc Abramowiczi I2000t IPen et al.l 120031: 
Bcgclman 2012; Mc Kinnev et al.ll2012D . in which case the 
total outflow rate integrated over disk radius is. 



Eo,, — 



2rfb V nn 



(3) 



where is the inner disk edge, either the black hole 
event horizon or the surface of the proto neutron star. 
For typical parameters we take rfb — 100 (~ 10®cm) 
and ^ = 0.1. These choices give an outflow energy 
Ew = eAffbC^ with e ~ lO^'^. For the parameters con- 
sidered here, other choices of s > give similar results. 
If the outflow is instead a jet launched from near the 
inner disk edge, the outflow energy is Ej = fiM^hC^, 
where a conventional choice is /3 ~ 0.1, although depend- 
ing on the accreted ma gnetic fleld geornetry t his value 
could be much larger ([McKinnev et al.[ 120121) . Since 
A/bh = A/fb(rtb/nn), the resulting energy injection. 
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Figure 1. Pressure (top), density (middle), and velocity (bottom) 
vs. radius at t = 100s, 300s, and 800s for an Eo = 1.2 X 10^^ erg 
explosion of a 25Mq zero metallicity progenitor using the fallback 
(lines) and piston (points) inner boundary conditions. The two 
methods are in good agreem ent away from the inner boundary, 
and the results are similar to IZhane et all ||2008| ') Fi Hure 1. The 
points have been down-sampled by a factor of 5 for clarity. 

eMfbC^ would be nearly identical to the case of a disk 
wind. This is the scenario discussed fo r failed supernova 
explosions bv iWooslev fc Heged (poTTI ) and QK12. The 
results presented below only depend on the energy injec- 
tion rate and thus are the same for either a disk or a jet. 
The results are also expected to be insensitive to whether 
the central object is a proto neutron star or black hole. 
In ijS.ll we discuss the dissipation of accretion energy in 
the infalling material and outgoing ejecta, and its impli- 
cations for the viability of wind and jet scenarios. 

We use e = 10~^ throughout, although we discuss disk 
formation and size in ij5.2l The outflow energy, either 
from a wind or jet, is then set by the fallback accretion 
rate. 

3. FALLBACK ACCRETION 

3.1. Numerical Hydrodynamics 

We estimate the fallback accretion rate by simulating 
supernova explosions using a ID lagrangian finite differ- 
ence hydrodynamics code. The code uses a st aggered 
mesh and artificial viscosity shock prescription (jCastoiH 
|2004[) . The artificial viscosity parameter is chosen to 
smooth shocks over ~ 7 zones. This is fairly diffusive, 
helping with code stability but gives nearly identical re- 
sults to much smaller coefficients. The Courant factor 
used is Aa;/csAt — 0.5, and the time step is set by 
the minimum required by any zone. The hydro code 
has been verified via comparisons to the Sedov- Taylor 
and ID shock tube problems, by verifying that pre- 
supernova stellar models with no explosions remain in 
hydrostatic equilibrium, and by comparing the solutions 




Figure 2. Numerical (lines) and semi-analytic (points) fallback 
accretion rates for three models from Table [T] For the small pro- 
genitor in S33, the asymptotic scaling M oc t~^/^ applies after the 
first ~ 100s. At much lower explosion energies (U45), the entire 
accretion rate curve is well described by freefall. In BSC progeni- 
tors (Z29), the accretion rate can be significantly enhanced at late 
times by material re-captured by the reverse shock formed at the 
H/He interface. 

for de-pressurized models with analytic freefall solutions 
(Eq. O. 

The initial conditions are taken from a wide 
range of pre- s upern ova progenitor star models from 
IWooslev et al.l ()2002f P1. Three sets of models are consid- 
ered: zero and solar metallicity progenitors with ZAMS 
masses of 11 — 4OM0, and \Q~'^Zq progenitors with 
ZAMS masses of 11 — 6OM0. The lower mass solar metal- 
licity progenitors retained their hydrogen envelopes and 
tended to be red supergiants {R ~ lO^^'cm), while those 
in the high mass range were bare helium or C/0 stars 
[R ^ lO^^cm). Low metallicity stars tended to be blue 
supergiants of smaller radii (i? ^ 10"'^^~^'^cm). Due to 
large uncertainties in prescriptions for semi-convection, 
convective overshoot, and mixing, we view the very low 
metallicity models as alternative outcomes for possible 
massive star progenitors rather than necessarily corre- 
sponding to extremely metal-poor environments. 

Explosions are simulated using a moving inner bound- 
ary ("piston," e.g. Wooslcy & Weaver 1995). For the 
first 0.45s the boundary moves inwards from the loca- 
tion where the specific entropy s = 4tor = 5x 10''cm, 
after which time it moves outwards at constant velocity. 
The inner boundary velocity is set to zero either after a 
specified amount of time or after the internal energy has 
changed by the desired amount. The resulting explosions 
are fairly insensitive to the piston velocity as long as it is 
large enough to deposit the desired amount of energy in 
a few seconds. We use the same number of radial zones 
for the hydrodynamics calculations as were used for the 
stellar evolution (~ 500— 1000), but verify that doubling 
the number of zones and interpolating using the near- 
est neighbor from the initial condition does not lead to 
significant changes in the evolution. 

An outflow boundary condition is employed by copy- 
ing the acceleration from the outer zone to a single ghost 
zone. The inner boundary condition can be either a hard 
(reflective) piston or inflow, depending on the time. Ini- 
tially, the piston is used to blow up the star, after which 
the inner boundary velocity is set to zero. To allow in- 

^ http:/ /homepages. spa. umn.edu/ alex/stellarevolution/data.shtml 
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Figure 3. Peak luminosity vs. time to peak for all events, measured from light curves calculated with the methods in Appendix A. 
The points are only shown when the peak luminosity is larger than the standard thermal supernova luminosity (Eq. |9]|- Luminous, long 
duration Type II events come from high energy explosions in massive BSGs, some with strong reverse shocks. Superluminous events with 
durations ~ 2 — 40 days come from weak explosions with low ejecta masses. Luminous Type I events are from weak explosions in compact 
initial stars. 



Table 1 

Sample Event Parameters 



Name 


MzAMS (Mq) 


MsN (M©) 


Z {Zq) 


Scxp (lO^'^ergs) 


Mej (Mq) 


Mrem (Mo) 


Wcj (km / s) 


Vj (km / s) 


Comparison 


S33 


33.0 


11.4 


1 


0.34 


2.1 


9.2 


2800 


24000 


SN 1998bw 


S39 


39.0 


8.49 


1 


0.21 


1.1 


7.4 


3000 


17000 


SN 2008D 


S390 


39.0 


8.49 


1 


0.055 


0.39 


8.1 


2600 


11000 


SN 2010X 


U45 


45 


44.7 


10-* 


1.0 X 10-^ 


0.45 


44 


340 


13000 


SN 2008es 


U60 


60 


59.2 


10-* 


1.0 


43 


17 


1100 


5500 




Z29 


29 


28.8 





1.1 


22 


6.7 


1600 


4700 






10" 10» 

Ejecta Velocity (cm / s) 



Figure 4. Peak luminosity vs. Vj for the events in Figure [3] 
The peak luminosity scales roughly with final ejecta velocity as 

flow, once the velocity of the inner zone drops below zero, 
it is copied to the inner boundary. This allows us to use 
the inner boundary to blow up the star and to record the 



accretion rate once material begins to fall back. When 
the inner zone passes through the radius corresponding 
to the assumed outer disk edge, rmin = 10* cm, its prop- 
erties are saved and it is removed from the calculation. 
The outside of the accreted zone then becomes the inner 
boundary for the subsequent evolution. 

We show in Figure [1] an example of the 1.2 x 10^^ erg 
explosion of a 25Mq, zero metallicity stellar progenitor. 
In this star, a significant density discontinuity at the he- 
lium/hydrogen interface (r ~ 2 x lO^'' cm) as well as the 
compact hydrogen envelope (p cx r"^'^) lead to a strong 
reverse shock forming at i ~ 20 s. The two-shock struc- 
ture can be seen clearly in the curves of p[r) and v(r). 
The filled circles show the results using a pure piston 
inner boundary condition, while the solid lines show the 
results for the "fallback" boundary condition (i.e., piston 
switched to inflow after the shock was initiated) . The two 
are in excellent agreement in the portion of the star with 
w > 0, but differ slightly in the inner regions, as pressure 
support slows the infall in the pure piston model. In the 
fallback calculation, the reverse shock turns around and 
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leads to a jump in the accretion rate at f ~ 2000 s. This 
calculation i s very similar to that shown in Figure 1 of 
iZhang et al.l (|2008f) . and the solutions are in qualitative 
agreement. The quantitative differences are likely due to 
a difference in progenitor models. 

We follow the explosions until late time {t = 10^ s), and 
calculate the accretion rate through the inner boundary 
from the properties of accreted zones. Sample accretion 
rate curves are shown in Figure [2] In some cases, we 
find large accretion rates (> 10~^ Mq s~^) for a week 
or so after the explosion. This late time accretion is due 
either to the fallback of stellar layers at large radii, or 
from the deceleration of inner layers by the reverse shock. 
The energy associated with accretion at these rates is 
sufficient to power luminous supernova light curves. 

3.2. Semi-Analytic Treatment 

The general behavior of the fallback accretion rate can 
be easily understood in the two limits where the ma- 
terial is either highly bound or mildly bound. For the 
highly bound material (i.e., those layers where the ve- 
locity following the shock propagation is much less than 
the escape velocity) the accretion rate can be estimated 
from the free-fall time, 



ts 



„3/2 



v/8GM(r) ' 



(4) 



from each radial and mass coordinate in the progenitor 
star: 



dM dr , , ^ dr 

M = — = 47rp(r)r^-— . 

dr dt ' dts 



(5) 



For an approximately power law density profile in a 
particular shell of a star, p{r) — Po{r /r^Y^"" , the en- 
closed mass is M{r) = po{ro) {r/ro)°' (for < a < 3), 
and the fallback accretion rate is: 



M = 



3 — a to 



3(a-l)/(3-a) 



(6) 



where <o = (2Gpo)"^/^ (cf. Eq. 2 of QK12). For a < 0, 
the enclosed mass is roughly constant, and the accretion 
rate is: 



M = 



8n porg 
3 to 



(2a-3)/3 



(7) 



where now to = 7rro^V\/2GM(ro). In this way, the 
freefall accretion rate is set by the density profile of the 
progenitor star. 

For the other limit of mildy bound material with Vcsc — 
V, the maximum radius, ri = ro(l — / Vosc)~^ i becomes 
much larger than the initial one, rg. Then the asymptotic 
fallba ck rate, M oct^ 5/3, applies (IMichellll988HChevalieri 
119891) . This asymptotic scaling applies at the latest times 
in all three curves in Figure [21 

Using the ballistics solution from lChevalieii (|1989[ ). we 
can bridge these two asymptotic limits to analytically es- 
timate the fallback accretion rate at all times for compar- 
ison with our numerical calculations. For each mass shell, 
the downstream shock velocity is taken from the analytic 



formulae in iMatzner fc McKei (|1999f) , which are typi- 
cally an excellent approximation to the numerical calcu- 
lations. Then the total fallback time for each mass ele- 
ment can be calculated from Eq. 3.7 of Chevalier (1989), 
and its time derivative is an approximate accretion rate. 
This assumes that pressure effects are negligible, which is 
incorrect. However, the true acceleration measured from 
the numerical calculations described below turns out to 
be roughly constant at half of the gravitational acceler- 
ation. 

This ballistic estimate reproduces the fallback accre- 
tion rate at all times in many progenitors. However, 
in some c ases (particular ly blue supergiants such as 
SN1987A. IChevahed 119891) the reverse shock formed at 
the hydrogen-helium interface is strong enough to decel- 
erate portions of the ejecta below the escape speed. This 
enhances the accretion rate at late time s, and c an sig- 
nificantly add to the remnant mass (Zh ang et al.l [2008). 
The reverse shock formation and evolution is analagous 
to that formed when the forward shock breaks out of 
the star and into the interstellar medium (e.g., iMcKed 
119741: lChevaliedri982[ ). As the simplest possible reverse 
shock prescription, we solve the strong shock jump con- 
ditions for the reverse shock velocity and the downstream 
velocity at the boundary of 100% helium and hydrogen 
layers: vrs — O-Gvq, where ?;o is the shock velocity. The 
reverse shock velocity evolves in time as the densities 
in both the expanding ejecta and unshocked hydrogen 
envelope change, and eventually it turns around. For 
simplicity, we ignore this and take vns to be constant 
at its initial value. Then the location of intersection be- 
tween ejecta and the reverse shock can be found, as well 
as the resulting ballistic t{M) for material that is re- 
captured after passing through the reverse shock. The 
reverse shock prescription is important for the Z29 curve 
in Figure [2j This approximate semi-analytic description 
does a reasonable job reproducing the numerical calcula- 
tions in all cases. The largest disagreement is in the re- 
verse shock cases, where the semi-analytic accretion rate 
overestimates (underestimates) the numerical results at 
early (late) times. For the remainder of the paper, we 
use the results from the numerical fallback calculations. 



4. POSSIBLE OUTCOMES 

We detail here the possible outcomes of supernova light 
curves powered by accretion energy. We first assume that 
a supernova explodes via the traditional core collapse 
mechanism, whatever that may be. For each progenitor, 
we ran explosions with a variety of energies, in the range 
10**^ to 10^^ ergs, in order to explore the full range of 
possible outcomes. Only explosions with positive total 
energy of non-accreted material at t = 10^ are consid- 
ered, and the resulting remnant vs. initial mass distri- 
buti on from thes e explosions is in excellent agreement 
with Zhang et all (|I008). 

The ejection of some stellar layers and the fallback of 
others is then calculated numerically as described in ^ 
which determines the energy input rate from fallback. 
We then calculate approximate one zone light curves us- 
ing the methods described in Appendix [XI For these cal- 
culations, we need the effective diffusion time through 
homologously expanding electa [Arnctt (1979.1 . 
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Figure 5. Comparison of fallback powered light curves (solid 
lines) from models U45 (purple), S33 (green), S39 (red), and S390 
(orange) with some observed supernovae. The parameters for these 
events are given in Table [T] The orange dashed curve assumes 
taff = 7 da ys. Data poin ts are taken from Gezari et al. ( 2003, 
SN 2008es). IMazzah et al.l ||200S . SN 1998bw and SN 2008D), and 
IKasliwal et al.l H20ia. SN 2010X1. 



3 Mk 
An vc 



3 (Mej+Mft)^ 



An 



VfC 



(8) 



where Mfb — ^ J Mfhdt is the total outflow mass, 
i?fb = eAffbC^ is the injected accretion energy, and 
Vf — yj {Esn + Efh)/ (Afcj + M[h) is the final ejecta ve- 
locity. Note that there is an ambiguity in determining 
Mfb, depending on the interpretation of the fudge factor 
^. If ^ indicates the fraction of outflow mass that inter- 
acts with the supernova ejecta, then the above expression 
for Mfb applies. If on the other hand, the mass transfer 
to the ejecta is more efficient while the specific energy of 
the outflow is lower, Affb could be significantly larger. 

We assume a constant opacity k — 0.2 g cm~^, appro- 
priate for electron scattering for fully ionized elements 
heavier than hydrogen. This is clearly a coarse approxi- 
mation, as the actual opacity will depend on the compo- 
sition and the presence of Doppler broadened lines. The 
effects of recombination on the opacity are, however, in- 
cluded in an approximate way (Appendix |^ . 

While our one zone light curve models account for the 
acceleration of the ejecta due to the input accretion en- 
ergy, they lack any information on the radial structure 
of the ejecta. The radiation hydrodynamical calculations 
of KBIO show that energy deposition at the base of the 
ejecta (in that case from a magnetar) blows a bubble in 
the inner regions, piling up material into a dense shell. 
We expect a similar effect in fallback powered SNe, which 
will likely also induce an asymmetry if the energy depo- 
sition is anisotropic. 

For each of the light curves, we measure the time to 
peak, tp, and the peak luminosity, Lp. The results are 
shown in Figure |3] for e — 10"'^. Each point represents 
a single explosion energy and progenitor model, color- 
coded by the radius of the pre-supernova star: red for 



m 10' 




100 200 300 400 
Days since explosion 



500 



Figure 6. Sample fallback powered light curves (top), and photo- 
spheric temperatures (middle) and velocities (bottom). The dashed 
curves in the top panel show the rate of energy injection from fall- 
back accretion. The temperature is estimated from the one zone 
model, while the velocity is taken to be the maximum of Vf and 
that calculated from the expanding ejecta. The temperature re- 
mains fixed at Tj during the plateau phase for events where hy- 
drogen is present. 

R > lO^^cni (RSGs), purple for lO^^cni < R < lO^cm, 
blue for lO^^cm < R < lO^^cm (BSGs), and green 
for R < lO^^cm (He or C/0 stars). This radius also 
corresponds to the zero age main sequence metallicity: 
solar without significant mass loss for RSGs, zero for 
BSGs, 1O~^Z0 for stars in between, and solar with large 
amounts of mass loss for compact He and C/0 stars. 
Events are only plotted if Lp is larger than the thermal 
supernova luminosity. 



Eq R 

td Vtd 



(9) 



The number of points is then set by the number of explo- 
sion energies and progenitor models, as well as the frac- 
tion of cases where that condition is met. The number 
of points does not represent an expected rate, since both 
the choices of explosion energies and progenitor models 
are arbitrary. 

Figure [3] illustrates the wide range of light curves that 
may result when fallback power is included. Many of the 
successful explosions with energies ~ 10^"'^ ergs lead to 
events with ~ 50 — 200 days, Lp ~ lO^^""*** ergss"^. 
The long durations are similar to those of Type II plateau 
SNe, and a result of the large ejecta masses and corre- 
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spondingly long diffusion times. The final velocities of 
these events are also fairly typical of core-collapse super- 
nova explosions (~ 3000 kms"^). This is because the 
amount of fallback is much less than the ejecta mass, 
so that fallback energy does not appreciably change the 
total kinetic energy of the explosion. For smaller ejecta 
masses, the fallback energy can dominate the total ex- 
plosion energy, significantly increasing the final velocity. 
The diffusion timescale therefore decreases with decreas- 
ing ejecta mass both from the smaller total mass and 
because of the increasing final ejecta velocity. 

These effects lead to a strong scaling of Lp with Vf, 
shown in Figure SI The roughly Lp cx vj dependence 
can be recovered by assuming the fallback energy always 
dominates the supernova energy (iJfb ~ vf), while the 
fallback mass contributes negligibly to the ejecta mass. 
Furthermore, the scaling assumes that the total fallback 
energy scales with peak luminosity {Lp ~ Ef^ ~ Wj), 
which is true if the accretion rate at late times scales 
with its integral over all times. The apparent maximum 
in ^/ ^ 3 X 10^cms~^ is from the case where the fall- 
back mass and energy dominate that of the supernova 
explosion: Vf ~ ^/ej^c ~ 0.01c for our standard pa- 
rameters. In the context of the simple outflow models 
described in Section [21 this maximum velocity scales as 
CX r^J . The considerable scatter in Figure g] is 
from the breakdown of the above assumptions. 

4.1. Candidate Events 

Different classes of progenitor stars lead to different 
outcomes in Figure [H First, solar metallicity RSG pro- 
genitors for the most part lead to relatively low luminos- 
ity events {Lp < 10"^ ergs s'^). At high ZAMS masses, 
these stars undergo substantial mass loss and become 
stripped He or C/0 stars. These progenitors can lead to 
events with tp ~ 20days, Lp ~ 10''^"'*^ ergs s"-^. These 
could potentially explain broad line Type Ibc GRB SNe: 
high velocities are a natural outcome of the injection of 
large amounts of fallback en ergy. Example fits a re shown 
in Figure E] for SN 1998bw f Galama et allllQQSl ) and SN 
2008D (Sodcrberg ct al. 200i). In the context of the col- 
lapsar model, this suggests that the central engine could 
be responsible for all of the observed properties: early 
time accretion leading to black hole formation, the GRB, 
and the initial supernova explosion; and late time ac- 
cretion powering the resulting light curve and the large 
expansion velocities. 

BSG progenitors lead to two classes of outcomes de- 
pending on the explosion energy. At low explosion ener- 
gies, they can produce luminosities Lp ~ io43-45 gj.gg g-i 
and peak times of ~ 2 — 40 days. The short du- 
rations are from the very small ejecta masses, Mcj ^ 
10~^ — IO^Mq, with nearly all of the star falling back. 
For ^ = 0.1 used here, the wind mass is comparable to 
the ejecta mass, and the injected fallback energy is much 
larger than the initial explosion energy. This leads to 
large final velocities and short diffusion times. Events 
with tp ^ 25 days can have light curve shapes very similar 
to observed luminous Type II-L events. A n example fit 
to th e superluminous Type II-L SN 2008es (|Gezari et al.l 
I2009D is shown in Figure [S] 

At high explosion energies, BSG progenitors lead to a 
range of long duration events with tp ~ 100 — 300 days. 



Lp ^ lO'*^"'*'* ergs s~^. The most luminous cases are ei- 
ther from very massive stars (> 4OM0) at low metallicity 
or from zero metallicity stars with strong reverse shocks. 
In both cases, the ejecta masses are ~ 10 — 4OM0 with 
low expansion velocities, w/ — 2000 — 6000kms~^. 

Subluminous Type I and II events are possible on a 
variety of timescales. As an example. Figure [S] shows 
a comp arison of a Type I e xplosion with the transient 
2010X (jKasliwal et al.ll2010[ ). The steep decay in this 
case requires that the accretion turn off about 7 days 
after explosion (see §5.ip . 

Approximate light curves from examples of each of 
these type of events are shown in Figure [5] along with 
photospheric temperatures and velocities. The model 
parameters are listed in Table [1] The photospheric tem- 
perature is taken from the one zone light curve calcu- 
lations (see Appendix |A| . The photospheric velocity is 
taken to be the maximum of f/ and the photospheric 
velocity in the expanding ejecta in the absence of in- 
jected accretion energy. For light curves with recombi- 
nation, the photospheric properties are meaningless after 
the plateau phase, since then formally the ejecta are com- 
pletely optically thin. For relatively short events, the ex- 
pansion velocities are high (> lO^cms"^, and the fallback 
energy sets the velocity since the ejecta mass is small 
(~ 1 — 2Mq). Much slower velocities occur in the longer 
duration events with large ejecta masses (~ 10 — 40Mq). 
The photospheric temperatures are very high at peak in 
ILL type events (~ 20000i4:). 

When recombination isn't important, the light curves 
are in excellent agreement with the semi-analytic for- 
mula in Eq. (jA6p for a power-law injection of energy 
with 71 = 5/3. This is because in nearly all cases the late 
time accretion rate falls as M cx t~^^^, while any energy 
injected on timescales < 1 day is lost to adiabatic expan- 
sion, so that its time-dependence does not influence the 
light curve. 

5. CAVEATS 

Gravitational energy liberated through fallback accre- 
tion at late times can power unusual supernova light 
curves (Figures [5] and [6|) . The calculations in this pa- 
per have made many simplifying assumptions; we dis- 
cuss here some of the uncertainties. We have treated 
the explosion of stars with crude ID hydrodynamic cal- 
culations using a piston. This method has frequently 
been used to simulate core coll apse supernova explosions 
and the resulting fallba c k Ce.g..lWooslev fc Weaverlll995l : 
iMacFadven et al.ll200H iZhang et al.ll2008[ ). and~the un- 
certainties in the numerically calculated fallback accre- 
tion rates are probably less than those in the outflow 
physics and/or parameters (e.g., e). The light curve cal- 
culations further assume simple one zone prescriptions 
for the bolonietric luminosity and photospheric temper- 
ature. More sophisticated techniques would be required 
for spectral calculations. 

The density structures of the pre-supernova models, 
which directly impact the fallback rate, depend sensi- 
tively on uncertain prescriptions for convection (semi- 
convection and overshoot) and co mpositional mixing i n 
stellar evolution calculations (e.g.. iWooslev et aI]|2002D . 
Probably a bigger issue is that the calculations here are 
based on a limited set of stellar progenitors, and ignore 
the effects of rotation and binarity, which may be very 
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Figure 7. Regions of 6j vs. t parameter space for model Z29 
where i) the outflow cannot escap e th e accreting material before 
depositing most of its energy (Eq. I12I I. ii) the outflow escape s th e 
outgoing supernova ejecta before losing most of its energy (Eq. I15I I. 
and iii) the energy deposited in the acc reting material by the out- 
flow exceeds its binding energy (Eq. I13I I. The remaining parameter 
space is where an outflow could plausibly power a supernova light 
curve without shutting off continuing accretion. For 9j < 10° , the 
outflow is arbitrarily changed from a wind {vj = 0.1c) to an ul- 
trarclativistic jet {vj ~ c). Constraints i) and ii) fix the range of 
allowed 9j for any disk formation time, ton, while constraint iii) 
sets the time at which fallback accretion will stop (toff). 

common in massive stars (e.g.. lSana et alJl2012t ). There 
may be additional variety in the range of possible fallback 
powered transients from stellar progenitors not consid- 
ered here. 

Further, we have assumed that the stellar material that 
falls back after the explosion has sufficient angular mo- 
mentum to form a disk, and that this disk can efficiently 
drive a massive wind and/or ultrarelativistic jet. These 
are both important open questions. The angular momen- 
tum distribution and surface rotation rates of massiv e 
stars remain highly uncertain (jWooslev fc Hege"ill2011[ ). 
Although previous studies have found prominent polar 
outflows from geometric a Uy thick black hole accretion 
flows^ (Stone et al.l 119991 : llgumenshchev fc Ab ramowic3 
I2000D . more recent calculations have found large-scale 
circula tions to be more common than unbound m assive 
winds ()McKinnev et al.ll2012HNaravan et"alll2012[ ). If so, 
ultrarelativistic jets may be a more natural explanation 
for injecting energy into the ejecta. 

Finally, we have assumed that this outflow can ther- 
malize in the outgoing supernova ejecta without expelling 
infalling material and halting accretion. We outline the 
requirements below to satisfy these assumptions, and es- 
timate in a few sample cases the required rotation rates 
for disk formation. 

5.1. Outflow Collimation and Interaction with Ejecta 

In order for a fallback accretion powered outflow (ei- 
ther a ultrarelativistic jet with speed Vj ~ c or massive 
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Figure 8. Sample fallback powered light curves (lines) for models 
(top panel) U45, S33, S390; and (bottom panel) U60 and Z29. The 
solid curves assume continued energy injection, while the dashed 
curves turn off at a range of times, t^ff = 0.3 — 4.0ttj. The dashed 
curves all assume constant opacity, and the strong effects of recom- 
bination can be seen on the light curves in the bottom panel. 
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Figure 9. Angular velocities required to form disks at radii at 
r = lO^cm (solid), lO'^cm (dotted), and 10^ cm (dashed) in a few 
models from material falling back at a range of turn on times. 
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wind with Vj ^ 0.1c) to power a supernova Hght curve, 
it must be able to: i) escape the remaining infalhng ma- 
terial ii) without unbinding it and iii) thermahze in the 
outgoing supernova ejecta. To order of magnitude, we 
assess the plausibility for this scenario as follows. 

Following QK12, we assume a magnetically-dominated 
outflow, collimated with an opening angle 9j. The prop- 
agation of the outflow through the remaining bound fall- 
back is similar to the propagati on of a jet through a host 
star during a long GRB fe.g.. iBegelman fc CiofEl IT989I : 
[Matzner 2003) . The speed of the head of the collimated 
outflow, Vh, is determined by pressure balance between 
the outflow and the host star (Eq. 4 of QK12): 



Vh 



1/2 



(10) 



where R is the maximum radius of bound material and 
M its total mass. The outflow escape timescale is then 
tcsc ^ R/vh- The outflow also drives a lateral shock into 
the surrounding material, whose speed is approximately 
(QK12), 



vi 



f 



l/4fll/2„l/4„3/4 



(11) 



where / is the efficiency of depositing outflow energy in 
the surrounding material. If the outflow is dominated by 
toroidal magnetic fleld (e.g., in a helical jet or outflow 
from a rotating disk), a typical value from num erical cal- 
culations is / ~ 0.03 ()Bucciantini et aLll2007D . For the 
outflow to escape, tesc should be shorter than the time 
for the lateral shock to envelope the star, ti ^ R/vi. For 
interactions with the remaining bound fallback material, 
this ratio is: 



0.6 



3 3 



1/8 



LjRih 



(12) 



In addition to requiring t; > tf,sc 
requires that the energy deposited. 



continued accretion 



/ 



ton 



Lj{t)dt, 



(13) 



should be less than the binding energy of the remain- 
ing fallback material, where ion (^off) is the time after 
explosion at which the outflow turns on (off). 

Conversely, for the outflow energy to be deposited ef- 
ficiently in the ejecta the outflow escape time should be 
shorter than the energy deposition time. Since the mech- 
anism for thermalizing the outflow energy and its associ- 
ated timescale are unknown, we can instead use the same 
comparison of ti and tesc as above. 

In this case, Vej is larger than Vh, and the escape time 
can be estimated from setting R = v^t and flnding when 



Vh 



(QK12): 



3Mejt;ej«jgf 

4L, 



Similarly, we can find the ratio U/tc 
assumptions. The result is: 



(14) 



under the same 



20 



m^3 



2/3 



(15) 



The requirements i) tcsc,FB/i;,FB < 1, ii) -Bdcp/i?grav < 
1, and iii) t;,cj/tGsc,ej < 1 amount to constraints on the 
outflow opening angle, 9j, and the time over which fall- 
back accretion can continue. Excluded regions of 6j vs. t 
parameter space from enforcing these constraints for the 
model Z29 are shown in Figure [Tj 

For the interaction of accretion energy with remaining 
fallback material, we find the maximum radius reached 
at time t by material that will ultimately accrete (i?fb), 
and its remaining total mass (Mfb) and binding energy 
(£^grav)- For the interaction of the accretion energy with 
the ejecta, we use M^] and Vf estimated at time t. 

The timescale constraints essentially place limits on 
6j for each type of outflow for all disk formation times, 
ion^ at small (large) opening angles, the outflow escapes 
(is captured). These ratios also depend on the other 
quantities, leading to differences between various mod- 
els. Generally, smaller Vf{t) leads to higher ejecta densi- 
ties and help to trap the outflow. Outflows escaping the 
ejecta before thermalizing could appear as long duration , 
high energy transients fQK12. IWooslev fc Hegeri[2011[) . 
Outflows trapped in the material still falling back would 
likely deposit energy there more effectively, either un- 
binding the material or prolonging its accretion to later 
times. 

At late times, the outflow will unbind any remaining 
material, shutting off further accretion. This is because 
the energy deposition into the accreting material at late 
times scales as Lt^sc t~^/^, while its binding energy 
scales as Mfb /rtb ^ t^^/"^ . Equating these gives the turn 
off time (toff) for each event. This turn off time tends to 
be shorter in higher energy explosions, since the bulk of 
the late time accretion comes from loosely bound mate- 
rial. Sample light curves for events where accretion shuts 
off are shown in Figure [8] for a range of ioff, assuming a 
constant opacity. Once the injected e nergy runs out, the 
light curve decays according to Eq. (jA5|) . but with an 
initial luminosity L(tos)- This may be particularly rele- 
vant for long duration transients in models like U60 and 
Z29, where the turn off time (~ 80 days for Z29) is likely 
to be comparable to the time to peak. 

Although these estimates demonstrate the plausibil- 
ity of accretion-driven outflows powering supernova light 
curves, detailed physical calculations will be required to 
assess this scenario accurately. Further, the statement 
that the outflow cannot escape the ejecta does not pro- 
vide an efficient means of thermalization, since we have 
assumed / = 0.03. Using a larger value of / ss 1 would 
shift the range of allowed opening angles to favor rela- 
tivistic jets, and lead to the outflow unbinding the ac- 
creting the material time at proportionally earlier times. 

The efficiency of thermalization depends on how ex- 
actly energy is transported from the accretion disk to 
the supernova ejecta. We have assumed that this mech- 
anism is a highly magnetized disk wind or an ultra- 
relativistic jet. If instead the wind is not highly mag- 
netized, a double (forward/reverse) shock structure will 
form when it catches up with the slowly moving inner 
layers of the supernova ejecta (KBIO). The situation is 
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analogous to the commonly case of supernovae interact- 
ing with circumsteUar material, only here the interaction 
happens inside, rather than outside the remnant. In ei- 
ther case, shocks should be efficient in thermahzing the 
kinetic ene rgy of the wind. Some rece nt semi-analytic 
(lBegelmaq,2012') an d numerical (McKi nnev et al.ll2012l : 
iNaravan et al.ll2012|) calculations of non-radiative accre- 
tion flows ha ve found large-scale circulations or convec- 
tive motions ([Lindner et al.|[2"011|) as well as or instead 
of outflows. This may also be a relevant mechanism for 
transporting accretion energy to large radius. 

If the accretion energy cannot efficiently thermalize, it 
will likely still lead to high ejecta velocities — Vf. In the 
case of an event like S33, this could still explain broad 
line Type Ibc supernovae: radioactivity would power the 
light curve and accretion energy would lead to the high 
observed photospheric velocities. This is also a possi- 
ble outcome of e arly time accretion onto a magnetar 
(IPiro fc Ottll201lD . 

5.2. Angular Momentum and Disk Formation 

Given the viable range of disk formation times for fall- 
back accretion powered supernovae 15.11 and the initial 
stellar radii accreting at those times, we can calculate the 
required angular velocity. Curves for models Z29, U45, 
and S33 are shown in Figure |9] for forming disks at radii 
from 10'^"^cm, or r ~ 1 - 100 for a lOMg black hole. 
Naively assuming rigid rotation, in all cases disks can 
form at the required times without exceeding breakup at 
the outer edge of the star. The required rotation rates 
essentially scale with explosion energy: for large explo- 
sion energies, the envelope is expelled, and larger rotation 
rates are required for the disk to form from material that 
was originally at smaller radius. 

Under this assumption, we can also calculate the max- 
imum disk size from fallback accretion, and the corre- 
sponding viscous time, tvisc ~ (R/ H)^a~^tdya, where 
H/R is the accretion flow scale height and a is the stan- 
dard dimcnsionlcss viscosity parameter in accretion the- 
ory (Shakura & Sunvacv 1973). Even with conservative 
assumptions [H/R = 0.1, a = 0.01), this timescale only 
becomes larger than the fallback timescale for the high- 
est rotation rates and early disk formation times < lO^s. 
This is because assuming rigid rotation, the total disk 
size never greatly exceeds its formation radius. 

Stellar cores spin up as they contract, and depend- 
ing on the efficienc y of angular m omentum transfer from 
magnetic torques ([Spruitl [2OO20 can transfer much of 
the core an gular momentum to the outer layers (e.g., 
IHeeer et al.ll2005l) . If this mass is retained, as in our 
models from low metallicity progenitors (e.g., Z29 and 
U45), it will likely form a disk upon fallback for modest 
ZAMS rotation rates. If instead this mass is lost (e.g., 
S33), insufficient angular momentum may remain to form 
a disk. If the red supergiant is in a binary system, tidal 
interactions may be an efficient means to spin up the star 
sufficiently to cause disk formation ev en if the envelope 
is los t during subsequent evolution (jWooslev fc HegeJ 
1201 It) . The latter scenario may be fairly co mmon, given 
the f requency of massive stars in binaries (jSana et al.l 
IMl) . These scenarios should be considered in more de- 
tail in future work. 

6. CONCLUSIONS 



The accretion power released when material falls back 
onto a compact remnant at late times could power un- 
usual supernova light curves. We have explored the con- 
sequences for a variety of progenitors and explosion en- 
ergies, using numerical calculations of the fallback ac- 
cretion rate and order of magnitude estimates of the 
resulting energy injection. While most of the fallback 
typically occurs at early times, it may be significant at 
late times in very massive stars, for low explosion ener- 
gies, or when a strong reverse shock forms at the hydro- 
gen/helium boundary. We have demonstrated that it is 
plausible that, under certain circumstances, the energy 
available from accretion could power an outflow which 
then thermalizes in the supernova ejecta. 

The events we have described are different and more 
diverse than what have previously be en studied as "fall- 
back supernovae". iFrver et al.l (|2009f). for example, con- 
sidered the case of massive star collapse in which most 
of the material fell into the central black hole and only a 
fraction was ejected. Because they also assumed that the 
surrounding medium was very dense and extended (due 
to mass loss prior to explosion) the supernova shock wave 
did not breakout of the circumsteUar gas until late times. 
The result was a dim, s hock-powered t ransie nt lasting 
from weeks to months. iMoriva et al.l (|2012f) similarly 
considered the case in which most of the star fell back 
and only a very small amount (~ O.IM©) was ejected. 
By assuming that this ejecta was enriched with ^^Ni, 
they found a brief and sub-luminous radioactively pow- 
ered transient similar to SN 2005E. Both of these pre- 
vious scenarios neglected the possible input of accretion 
energy from fallback (i.e., they assumed e = 0). As we 
have shown, accretion may re-energize the ejecta at late 
times and hence power much brighter emission. 

The power from fallback accretion may be relevant 
for explaining recently discovered classes of peculiar su- 
pernovae. These may include the Type IIL supernovae 
that are extr emely luminous and of relatively short du- 
ration (e.g., iGezari et all 120091: iMiller et all 120091 ) as 
well as those that are moderately bright and of very 
long duration (e.g. [Miller et al.l [20Tol: iRest et "all 120111: 
iChatzopoulos et aLl 120111 ). Several of the observed 
Type II events, however, also show narrow hydrogen 
emission features in their spectra, indicating that inter- 
action with a dense circumsteUar medium is occurring 
and may be responsible for the luminosity. 

Many of the predicted Type II events with ^ 10^^ erg 
explosion energies have very long times to peak (100—200 
days). Accretion energy is likely to unbind the remaining 
infalling material on a comparable timescale, turning off 
the power source for the light curve ( §5.1|) . We therefore 
predict that these events could be seen as very bright 
Type II supernovae that disappear suddenly from view. 
In general, late time turn off would be an observational 
signature of fallback accretion powered supernovae. 

Fallback accretion could also power very bright Type I 
events. The models considered in this paper reached 
peak luminosities of ^ 10'*'^ ergs s^^, similar to the 
broad-lined SNe Ic like SN 1998bw. If the accretion 
efficiency is assumed to be higher than our fiducial 
case, it is possible for some events to reach luminosities 
> 10**^ ergs s~^, in which case fallback could power the 
super-luminous hydrogen poor events such as SN 2005ap 
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(iQuimbv et al.|[20Tll) . 

Another effect that may produce super-luminous 
events hke SN 2005ap involves mass loss. Some events 
considered here (e.g., Z29) are brightened considerably 
by enhanced accretion from material decelerated by a 
reverse shock forming at the H/He interface. A similar 
outcome could occur in both Type I/II events where the 
progenitor has experienced considerable mass loss shortly 
before explosion. In this case, the reverse shock would 
be formed when the outgoing shock wave reaches the 
interface between the progenitor star and the massive 
wind or ejected shell. The subsequent inward propaga- 
tion of the reverse shock could lead to order of mag- 
nitude increases in the fallback accretion rate at later 
times. Interaction of supernova ejecta with circumstellar 
shells at radii ~ 10^^ cm is commonly considered to ex- 
plain superluminous supcrnovac via thermalization of the 
kinetic energy (e.g.. iGal-Y am 2012,). Surprisingly, inter- 
action with circumstellar material at much smaller radii 
(~ 10^^ — 10^^ cm) may also lead to an super-luminous 
event, but by very different means - by enhancing fall- 
back accretion and feeding the central compact object at 
late times. 

Also of potential relevance to the fallback scenario are 
dimmer supernovae that decline very rapidly after peak 
(e.g., SN 2002bj an d SN 2010X. IPoinanski et al.l[2010l: 
iKasliwal et al.|[2010D . The short duration of these events 
makes it difficult to explain them as radioactively pow- 
ered transients. In the case of SN 2002bj at least, the 
mass of ^^Ni inferred from the light curve peak exceeds 
the total ejecta mass inferred from the light curve dura- 
tion (diffusion time), seemingly ruling out a radioactively 
powered events. In the fallback scenario, short lived tran- 
sients are possible, especially if the energy injection from 
accretion cuts off the fallback abruptly (Figure [8]). 

The true range of possible light curves powered by 
fallback accretion is likely much larger than shown in 
Figure [3] Those events are limited in both the variety 
of progenitor models and by our neglect of the fallback 
accretion physics. The latter could conceivably lead to 
variations in e in either direction: smaller disks and/or 
more efficient thermalization of the outflow could lead 
to higher peak luminosities > lO''^ ergs s~^. Conversely, 
lower efficiencies could help explain a wide r variety of 
sub-l uminous supernovae (e.g., SN 2008ha, iFolev et al.l 

Further modeling is needed to identify the observa- 
tional signatures of fallback accretion powered super- 
novae, and determine how we might distinguish these 
events from other means of generating unusual light 
curves. Possible signatures include a tail with L cx t^^/'^ 
at late times compared to the diffusion time (but be- 
fore the ejecta become completely optically thin), some- 
what different from those from radioactivity or magne- 
tar spindown. More promisingly, at late times (^ 100 
days) it seems likely that the accretion energy will un- 
bind the infalling material (§ 15. ip . shutting off accre- 
tion and leading to a sudden decrease in luminosity. If 
instead accretion continues for decades after the explo- 
sion, th e black hole could em erge as an observable X-ray 
source (iBalberg et al.l 120001) . as has been suggested for 
SN 1979C (Patnau deet al.ll2011[ ). 

The conditions for fallback to influence the supernova 



light curve are apparently quite special, as the scenario 
requires sufficient angular momentum to form a disk and 
an evolution that permits fallback to persist long enough 
to drive energetic outflows at late times. Such a conflu- 
ence of factors may be rare in the Universe. On the other 
hand, observational surveys show that the rate of pecu- 
liar SNe - in particular the rate of very luminous ones 
- is only a small fraction that of standard core collapse 
events. It is possible that fallback power plays a role in 
some of these spectacular events. 
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APPENDIX 
LIGHT CURVE MODELING 

KB 10 described a one zone diffusion estimate for bolometric supernova light cur ves powered by an inj ection of energy 
with arbitrary time-dependence, H(t). The argument follows along the lines of I ArnettI (|1979l . [r982[i . As the ejecta 
expand, energy is lost both due to escaping radiation (L) and adiabatic losses from expansion: 



dV 



-P- 



H - L 



at ^ dt 

Assuming optical depth r 3> 1, the diffusion equation gives an approximate relationship between Ei^t and L: 

^ _ L _ c Ei,-,t/V 
47ri?2 ^ 3kp R 

With the definition of the diffusion time (Eq. [8|) and assuming V oc t^, Eq. (jAl|) can be re- written as: 



d 
di 



The general solution for L{0) = is: 



L{t) = ^H{t). 



dtte^^/^^^'^ H{t). 



For an initial shock energy H = Eod{t — to), this gives: 



Lit) 



td td 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



in agreement with Eq. (|9|) with to = R/vgh modulo the exponential factor of order unity at peak. 

For accretion powered light curves, a power law form, H — _Lo(i/io)~" with n = 5/3 provides an excellent ap- 
proximate description for the numerical light curves integrated with Eq. (jA4|) . The semi-analytic solution for t > to 
is: 



L{t) = Lo 



-eiiti 



V2 



7 1 



n 

2' 2tl 



7 1 



n 
2' 



t^ 



(A6) 



where "f{s,x) is the lower incomplete Gamma function. The incomplete Gamma function is complex for negative 
arguments. Since the observed light curve and the integral in Eq. (jA4|) are real, the imaginary part in Eq. ()A6[) 
vanishes. In the special case of constant energy injection {n ~ 0), the solution is (cf. Eq. 13 of KBIO): 



L{t) = Lo 



1 _ e-^t^-tl)mi 



(A7) 



for t < tos, where iJ = for i > tos- 

This light curve estimate assumes a constant opacity. A different limit occurs when the outer portion of the ejecta 
drops below the ionization temperature and recombines. The opacity drops suddenly in the recombined material, and 
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the effect is that of a recombination wa ve passing through the ejecta. This effect significantly alters the light curve 
evolution of Type II-P supernovae (e.g., lPopovl[r993t iKasen fc Wooslevl[2009l) . 

During the passage of the recombination wave through the ejecta, the photosphere remains at the ionization tem- 
perature, Tf. The luminosity can then be calculated from the time-dependent photospheric radius: 

L = AT:Rl{t)aTf, (A8) 

where Rp(t) ~ Xi{t)vt and Xi is the dimensionless position of the photosphere in the expanding ejecta. We can write 
the equivalent of Equation (|Al|) for the evolution of the internal energy of the ionized region, eVi — A/'ine{xiVt)'^: 

Ide A dx^ A H-L 

Idt^x^^^l^^vT' ^ ' 

We again use the diffusion equation to write L in terms of the internal energy, except now using Rp instead of R. 
Finally, we equate the photospheric luminosity with that from diffusion in the ionized region, which gives an expression 
for e in terms of Xi . The result is a non-linear first order differential equation for Xi (t) : 

dt 5t 5tjx, 5xlt \A^vHluTf ) ■ ^ ' 

In the absence of heating {H = 0), Eq. (jAlOp is similar to Eq. 14 of iPopovl ()1993[ ). except with slightly different 
numerical coefficients. In this case, the analytic solution for the luminosity starting at time tt, such that Xi(ti) — 1, is: 



L{t) = ATiaTfv^ 



7tl) 7t^ 



(All) 



In general, we calculate the luminosity assuming constant opacity using Eq. (jA4l) . Then, the approximate one 
zone photospheric temperature is given by aT^ = L/AttvH'^. When this drops below T/, we numerically integrate Eq. 
(jA10[) for Xi{t), and then calculate L{t) = Airv'jt'^ xf aTf . The recombination wave can significantly increase the peak 

luminosity in hydrogen rich progenitors (see Section 14]). More accurate radiative transfer calculations would likely find 
smoother light curves than those estimated from this one zone approach. 



